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Abstract 

In this paper we investigate complex solutions of the Bethe equations in 
the two-particle sector both for arbitrary finite number of sites and for the 
thermodynamic limit . We find the number of complex solutions (strings) 
and compare it with the string conjecture prediction. Some simple properties 
of these solutions like position in the spectrum, crossing of levels, connection 
to the ground state and transformation to the real solutions are discussed. 
Counting both real and complex solutions we find expected number of highest 
weight Bethe states. 

I. INTRODUCTION 

Among integrable spin chains those invariant on quantum groups have attracted consid- 
erable interest. The simplest among these are SU(2) q invariant chains. Open chains have 
been considered for spin one half spin one and higher spin [|[. Generalisations to 
other groups have been also investigated ||], |5|]. Closed chains have been introduced later 
because parallel requirement of quantum group invariance and generalised translational in- 
variance required introduction of a nonlocal term in the Hamiltonian ||. These chains were 
shown to have interesting properties. SU(2) q invariant closed chain has ground state with 
vanishing or nonvanishing spin depending on the value of the coupling constant M. Central 
charge was found and it was shown that in particular points of the coupling constant it cor- 
responds to central charge of minimal unitary series J7|, ||. Its excited states and operator 
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content was also found 0. Recently, it was argued that this model is related to interesting 
field theories [T(J and in particular to the Liouville theory with imaginary coupling [pTTfl . A 



common approach to investigation of integrable models in general and of this model in par- 
ticular is the Bethe ansatz method ||12|| , fll3| . This method leads to a set of transcendental 



equations for momenta of quasiparticles. One class of the solutions has real quasimomenta 
that can be found by numerical iteration. However, there are also solutions with complex 
quasimomenta. Searching for the latter class of solutions one usually makes the so-called 



string conjecture |14[], [15|. In particular, it has been used in investigations of some quantum 
invariant models [|T^], [17] and also in a recent investigation motivated by the relation of 
such models to Liouville field theories |18[]. However, it was recently pointed in the context 
of XXX and XXZ models that string conjecture has exceptions fi~9H25l1. That means those 
solutions to Bethe Ansatz equations are not yet completely understood. For that reason in 
this paper we investigate complex solutions for the closed quantum invariant spin chain. In 
this investigation we shall not use the string conjecture. 

For simplicity, we shall investigate the sector where spin is lower by two units than the 
maximal spin. We shall find the number of complex solutions for arbitrary number of sites 
N. In the thermodynamic limit, this number is found to be of the order N, as expected on 
the basis of the string conjecture. Next correction will differ for a finite number from the 
string conjecture prediction. If the coupling constant tends to the Bethe point (cos</> — ► 1) 
in a sufficiently fast manner, than we shall have an infinite number of exceptions to the 
number predicted by the string conjecture. The model in this point coincides with the 
XXX model and this is consistent with results of [^], | 2"3fl . Further, we shall find that the 



energy distribution of complex solutions show simple features and that they are located in 
a narrow energy band. In particular they are located on the top of the spectrum near the 
antiferromagnetic point and on the bottom near the ferromagnetic point. Energy levels cross 
each other as coupling constant changes. There are no such crossings in the XXZ model. 
We also find that solutions of the one of the two classes of bound states (strings) evolve 
in real solutions in the special points for the coupling constant where we know that the 
representation theory is not isomorphic to usual SU(2). It will turn also that at least one 
of overall ground states evolves in the string state. This happens in a region of coupling 
constant where it is not more ground state. 

As already anticipated we shall consider the Hamiltonian 

JV-l 

H = Nq-Y,Ri-Ro (1-1) 

i=l 

Rq = GRn^iG 1 (1-2) 

G = R1R2 — Rn— 1 (1-3) 
where Ri are 4x4 matrices given with 

Ri = a+*. +1 + oral, + (ajVf +1 + l) - (<t? - af +1 - 2) . (1.4) 

Here q denotes a parameter which lies on unit circle 
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q = e 



i<j> 



;i.s) 



and cf) will be called coupling constant. The operators Ri satisfy Hecke algebra 



(q - q- 1 ) R t + 1 



;i.6) 



As a consequence 



RiRi+lRi — Ri+lRiRi+l- 



[G, H] = 0. 



;i.t) 



;i.8) 



In addition it was shown that this Hamiltonian is invariant on the SU(2) q symmetry group 
whose generators are given with 



5 3 = -Vl 

2 frf 

i=i 



07 



;i.9) 



-a 3 /2 



y a 3 /2 



i=l 



\s + ,s- 



q 2Si - q 253 
q - q- 1 



;i.io) 



;i.ii) 



Hamiltonian of the model is highly nonlocal but due to Hecke algebra it is still integrable. 
It can be diagonalized e.g. with the coordinate Bethe Ansatz method 0]. The energy 
eigenfunctions of the spin s 
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with M spins down can be written as 



\^m) = E tpM(ni,...n M ) \ni...n M ). 

Kni<n2...<njif<)V 



The ipM functions are given with 



;i.i2) 



1.13) 



ip M {ni, . ..n M ) = J^exp 
p 



M 



E^% + o E *Pjfi 



i<j'<'<a/ 



;i.i4) 



Here the sum runs over the elements of the permutation group Sm- The phase shifts $jj 
have the following simple expression 



The quasimomenta ki, i = 1 . . . M form a solution of the Bethe Ansatz equations 

M 

Nki + (2M — N — 2) + — 27rAj i = l...M. (1.16) 

i=i 

The M Bethe numbers Ij, % = 1 . . . M are half integers (integers) for M even (odd). In terms 
of quasimomenta ki,i = 1 . . . M the energy E and generalized momentum P read 

E = 2^(cos0- cos^), (1.17) 



AY 



P = ^fc i -0(iV-M-l). (1.18) 

i=i 

Operator G is then 

G = e- iP . (1.19) 

Due to SU(2) q symmetry one can use its representation theory to classify the states. In fact 
for generic q we have the same multiplet structure as for usual 577(2). However, for 

q v = ±1 p integer (1.20) 

additional degeneracy occur [l|], ||. In particular representations of spins j = j + np and 
j = p — 1 — j — np mix. Here n is an integer. In order to get representations isomorphic to 
SU{2) one has to exclude also 

j = np -1/2. (1.21) 
Only the remaining representations (called 'good' representations) with 

j < V -^~ (1-22) 
are isomorphic to 577(2). The parameter q which is a root of unity can be written as 

q = e l ~. (1.23) 



It was shown [|] that for p' = 1 the representations are unitary. We shall find afterwards 
that the points defined with ( |1.23| ) will play a role in the evolution of complex solutions into 
real ones when we vary the coupling constant. Let us for the moment concentrate to the 
generic sector of the parameter q. Due to the already mentioned property that in this case 
we have the usual 577(2) multiplet structure it is sufficient to find the highest weight states. 
Other states can be constructed J7| by application of the generator J~. It is known that the 
highest weight states correspond to sets of {k 1 }, i = 1 . . . M where the quasimomenta satisfy 

k + (p. (1.24) 
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In fact we can identify in advance Bethe numbers which lead to the non-highest weight 
solutions. Let us assume that for some % = I we have 



k = 0. (1.25) 
A straightforward calculation using ( 1.15Q shows that 

{k = 0, kj) = 7T - 20. (1.26) 
From Bethe equation Ql.16 ) for % = I we obtain 

M-l 

A* = — (1.27) 

The highest weight solutions will be obtained by excluding ( 1.27j ) from the choice of the 
Bethe numbers. The ground state for the Hamiltonian ( |1 . 1| ) was found for the whole interval 
< < it. In fact, contrary to the the XXZ model, the spin zero state is the ground state 
only in the | < < it region. In the rest of the interval there is a subregion for each spin 
where the ground state has just that particular spin. More precisely, the total spin s of the 
ground state depends on the coupling constant according to 

J = for -<0<7T 

= * for 2(£lj*^S (L28) 

N 7T 

J = — for < < — . 
2 ~ N 

The Bethe numbers which give ground state are given with 

/ M-l\ /M-l\ 

A =(-_)-!,...,(_ )-l. (1.29) 

and 

M = — -s. (1.30) 

It also turns out that the corresponding momenta are real. The natural question arises 
which role play complex solutions of Bethe equations. 



II. COMPLEX SOLUTIONS IN TWO PARTICLE SECTOR 



In this section we shall investigate complex solutions of the Bethe equations ( |1.16| ) with- 
out assuming the string conjecture fl4||, [15]. For simplicity we shall work in the M = 2 



sector. In that case the Bethe equations ( |1.16|) take the form 



Nki + $i j2 + 0(2 — N) = 2ttAi, 
Nk 2 - $1,2 + 0(2 - N) = 2tt\ 2 . 



(2.1) 
(2.2) 
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Here we want in particular to look for complex solutions. Due to reality of energy and 
generalized momentum, k\ and k 2 have to be complex conjugates of each other 



k\ — kf ~\~ ik'i^ 

We can express k r and fcj by taking sum and difference of equations 

1 



N 



[tt(\i + A 2 ) - 0(2 - AT)] 



1 and (12 



(2.3) 
(2.4) 



(2.5) 



iNki = 7r(Ai — A2) — 2 arctan 



cos0sin(iL- 



cos k r — cos cos(iki) 



With the help of the identity 



2i 1 — iz 



and exponentiation of (|2.6|) one can obtain 



sinh(fcj(y — 1)) cos/c r 



sinh(/cjY 



cos < 



Ai + A2 odd, 



(2.6) 



(2.7) 



(2.8) 



cosh(/cj(y — 1)) cosfc r 



cosh(fc. 



2 - 



COS( 



Ai + A 2 even. 



(2.9) 



In the further analysis we shall call the solutions of equation ( |2.8| ) s-solutions and the so- 
lutions of equation ( |2.9| ) c-solutions. The left hand sides of both equations (|2l 



if) are 



monotonously decreasing functions so we shall have a solution for ki for any k r for which 
cos k r is in the interval 



if cos > (0 < (j) < f 



s-strings 
c-strings 

and 



< cosAv < cos 0(1 — — ), 
< cosAv < cos 6 



(2.10) 
(2.11) 



s-strings : cos0(l — — ) < cos£; r < 0, 
c-strings : cos < cos k r < 



(2.12) 
(2.13) 



if cos0 < (| < (p < 7r). Further, as long as ki 7^ 0, k\^ 7^ <f) an d thus this solution 
represents the highest weight state. Now we can proceed to find the number of complex 
solutions. As a first step we shall consider leading order in A^ when the inequalities (|2.10 ) 
and ( |2.12| ) are identical. In this case, for a given coupling constant 0, the number of complex 
solutions will depend on the number of values that k r can take. In order to have one to one 
correspondence between k r and COS tllG sum of Bethe numbers can take 2A^ — 1 different 
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equidistant values . The interval in k r for which there are complex solutions is 2 (~ — 0^ . 
As a result the number of the string solutions in leading order in N is 

J_(2iV- 1)^-20). (2.14) 
Ztt 

Maybe we shall remark that if we insist to have k r between 7r and —it we should not take 
Bethe numbers symmetrically around zero due to the term proportional to in equation 
( p.5[ ). However, this does not affect the counting argument. To determine the number of 
complex solutions more precisely, we have to take into account subleading orders in N in 
relations ( |2.1U| ) and ( [2.12 ). The allowed interval for the real parts of s-solutions is smaller 



than the corresponding interval for c-solutions and there are no s-solutions for an interval 
in k r of length 

COS 

Sk r = 4arcsin -, , (2-15) 

„_ . / arccos cos rf>+arccos( cos Ml — -Fr) ) x 

N sin ' * — - — ElJ - 



The number of solutions, which were overcounted in naive formula ( |2.14|) , is integer part of 



2N-1 . cos v 
n = arcsm -, . . (2-lo) 

7f / arccos cos <j!>+arccos( cos <j>{l — 4t) I 

JM sm 1 £L - L 



This correction is a finite number even in the thermodynamic limit (N — > oo). However, if 
| cos 0| = 1 correction is infinite and goes as y/N. This is true if | cos0| — > 1 sufficiently fast, 
more precisely if 

iV(l - | cos0|) a = const, and a > 1. (2.17) 

This is consistent with the fact that = or = pi corresponds to XXX chain for which is 
known that this correction is infinite. This correction represents at the same time a violation 



of the string conjecture which was found previously in XXX p2| , |p3|| and XXZ chain |24 
f25| . Disappearance of s-solutions (strings) is followed by appearance of real solution with 
two close quasimomenta. These quasimomenta can be found near the number of sites where 
s-string disappeared by numerical iteration of ( |2.1| ) and ( |2.2p with the same Bethe numbers 
Ai and A2. Real solutions of Bethe equations with two identical Bethe numbers also represent 
violation of the string conjecture. Evolution of s-strings into real solutions with identical 
Bethe numbers can be followed for fixed coupling constant by increasing the number of sites 
N. Exapmle is given on Fig.l for = 0.32. 

III. PROPERTIES OF COMPLEX SOLUTIONS 

Next interesting question we would like to ask is how are complex solutions (bound 
states) distributed on the energy scale. From the relations (|2.9p, ( |2.8| ) and ( |1.17[ ) we see 
that complex solutions are confined in the energy band 

< ^(c-strings) < 2cos0, (3.1) 

8 cos „ . . 
— y < ^(s-strings) < 2 cos0. (3.2) 
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The left sides of the inequalities correspond to the points (in coupling constant 0) where 
two complex quasimomenta collide on the real axis and the complex solution becomes real 
solution (decay of the bound state). The right sides of the inequalities correspond to the 
points where ki tends to infinity and so the localisation of two overturned spins in bound 
state tends to infinity. This can be seen from the form of the Bethe wave function ( |1.14[ ). 



The energy band ( |3.1|) is generally narrow compared to the overall spread of energy when 
all solutions are included. As an illustration of M = 2 spectrum, figures are given for 
N = 6,N = 8 and N = 10 (Fig.2, Fig.3, Fig.4). We see that all bound states (strings) 
disappear near the 'free theory' point = |. Energies of the string solutions are on the 
top of the spectrum near the Bethe antiferromagnetic point = 7r and on the bottom of 
the spectrum near the Bethe ferromagnetic point = 0. There is an interesting question 
connected to the nature of the overall ground state. We know that for any finite and even 
N the total spin s of the ground state depends on the value of the coupling 0. In particular, 
ground state has spin zero for ~ < < ix and spin s for 2 u+x) — ^ — Is" ^ ne quasimomenta 
are known and they are all real. We find that beyond these intervals there is a region of 
values where these solutions become complex. To see this feature we remind that due to the 
relation M = y — s and fll.2S| ) the overall ground state is in the M = 2 sector for 

^ <0<-T7^- T . (3.3) 



N-2 ~ T ~ JV-4 

In particular for iV = 6,8,10 these intervals in coupling will be [f ,§],[§, f] and [|,f] , 
respectively. Following this state we see that at certain points outside the above intervals 
the quasimomenta become complex. From the choice of Bethe numbers ( 1.2 9|) we see that 



this state will become complex c-solution and the transformation of this real state into string 
will happen in the E = point. Generally, we can find points in coupling constant where 
the complex solutions will become real. In the transition points 

cos0 = cos/c r c-strings (3.4) 



N 

For c-strings this is satisfied for 



cos0 ( 1 — — ) = cos&v s-strings. (3.5) 



" -[2,...,iV-4,iV-2] O<0<- (3.6) 



r 2(iV-l) L '"'' ' J ~ r ~ 2 

0=— — — r[JV,...,2iV-6,2jV-4] -<0<tt. (3.7) 
r 2(N- 1) L ' ' ' J 2 ~ y ~ K ' 

Now we see that the state, which is the overall ground state in interval ( |3.3| ), after a change of 
(n-i)(n-2) m ^ becomes c-string. It is interesting that the points ( |3.6|) and ( |3.7|) correspond 
to the points ( |1.23|) where the the representation theory is no more isomorphic to SU(2) 
and where several multiplets merge together forming indecomposable combinations jlj . The 
points where s-strings disappear are not of this form because of the correction factor ( 1 — 
in (U). 
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Crossing of energy levels of complex solutions with the change of coupling constant is one 
of the features of this model that make it different from the XXZ model. This is illustrated 
on Fig. 5 and Fig. 6 on which we follow energy levels of M = 2 strings for both models and 
for number of sites iV = 15. The number of strings is the same in both models. However, in 
XXZ model there are complex solutions with ±k r which are degenerated. Presence of a term 
linear in <p in equations (|1.16|) , which can be interpreted as a coupling constant dependent 
toroidal twist in XXZ model, removes this degeneracy and causes crossing of energy levels. 

IV. REAL QUASIMOMENTA AND COMPLETENESS OF SOLUTIONS TO 

BETHE EQUATIONS 



Now we want to enumerate all real solutions of Bethe equations in two-particle sector. 
We start again from equations (|2.1|) and (2.2). After manipulating their difference and sum 
we obtain the following equations for k = k± + k 2 and k\ — k 2 



k 



N 



br(Ai + A 2 ) - 0(2 - iV)] 



(4.1) 



Sill 1 



(fci-fa) (N 



1)) 



sin 



(fa-fcj) AT' 



cos • 



COS ( 



Ai + A 2 odd 



(4.2) 



cos(^l(f -1)) 



COS! 



(fci-fc 2 ) N- 



cos ■ 



COS i 



Ai + A 2 even. 



(4.3) 



We have to notice that not all 2A" — 1 different values of / = Ai + A 2 give different 
solutions. Changing / by AMs equivalent to changing one quasimomentum by 2tt which 
results in change of sign of right hand sides of fl4.2Q and ( fl.3| ) . This reduces the number of 
possible values of / to N, which in turn can be chosen to give positive values of right hand 
sides of ( [4.2[ ) and (p~3|). The left hand sides of the equations ( |4.2j ) and ( |4.3| ) are periodic 
functions. Thus in principle, for each of A" different fixed values of the right hand side 
one can count number of solutions by counting number of intersections. These periodic 
functions are given on Fig. 7 for A^ = 9 and A" = 10. For A" odd we find that the number of 
intersections is 



N 



'cos I 

COS (j) 



1 + 



A^-3 




I even 



(4.4) 



'cos I 
cos (b 




I odd 



(4.5) 



while for A^ even it is 

Nt = 




1 + 



N — e d cos| 



I even 



COS( 



(4.6) 
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ATO N -2 /cos | i 2\ A^ — 4 / 2 cos§\ T 

N° = ——9 2 - - 1 + — + — — 9 1 - — 1 I odd. 4.7 

2 \cos0 TV/ 2 \ N cos0/ v 7 

The number of complex solutions is just equal to the second 9 function in the expressions 
above. If we take into account that for N even there are y even and y odd values of / and 
for iV odd ^-=^ even values of I and odd values of /, we find that for any number of 
sites N and coupling constant <f> there are 

N(N — 1) , , 

1 2 ^ (4-8) 

solutions to the Bethe equations in the two-particle sector. Among these solutions there are 
A" previously identified non-highest weight states, which have one of quasimomenta equal 
to 0. Finally, we obtain 




highest weight Bethe states. 



V. CONCLUSION 

In this paper we have investigated complex solutions of the Bethe equations for the M = 2 
sector of the SU(2) q invariant closed spin chain for arbitrary number of sites and coupling 
constant. We find that some properties of these solutions are similar to the properties of 
complex solutions for the normal XXZ chain and some of them are not. In particular we find 
the number of complex solutions both for finite A" and in thermodynamic limit. This number 
differs from the number predicted by the string conjecture for a finite number of solutions. 
However, if <fi — > sufficiently fast compared to the number of exceptions becomes 
infinite and goes as \fN . These properties are essentially the same as in the case of XXZ 
chain. One can follow evolution of string solutions and their disappearance with decreasing 
'coupling strength' |cos0|. One class of strings (c-strings) disappears in points where the 
representation of the SU(2) q is no more isomorphic to SU(2). By increasing number of 
sites for fixed coupling constant, strings of the other class turn into real quasimomenta that 
can be found by iterating Bethe equations with two identical Bethe numbers. This is again 
violation of the string conjecture. The overall ground state, which is always real, is in the 
M = 2 sector for coupling constant — — ]v~br ^ n ^ s state becomes complex solution 
for <fi < j^zi- The energy dependence of string solutions shows some simple features. The 
strings are found in a narrow energy band and are located on the top of the spectrum near 
the antiferro magnetic point and on the bottom near the ferromagnetic point. Their energy 
levels cross each other with the change of coupling constant, which is not the case for the 
XXZ chain. Finally, we find the number of real solutions. The number of all solutions of the 
Bethe equations is [^j . Among these we identify A^ non-highest weight states. That leads 

to the [Z) — y[ \ highest weight Bethe states. 
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FIGURES 



FIG. 1. This figure shows dependence of real part of complex s-solutions on numbers of sites 
N for <fi = 0.32. It clearly illustrates transmutation of one complex solution in real solution (two 
quasimomenta) for a given critical N. These two real quasimomenta correspond to same Bethe 
numbers and are obtained by numerical iteration of equations Q2.1| ) and (|2.2j). 

FIG. 2. Spectrum of M = 2 sector for N = 6. Energies of both real and complex solutions of 
the Bethe equations are plotted. 

FIG. 3. Same as Fig. 2 but for N = 8. 

FIG. 4. Same as Fig. 2 but for N = 10. 

FIG. 5. Dependence of energies of complex solutions for SU (2) q invariant spin chain and num- 
ber of sites N = 15 on coupling constant (j>. 

FIG. 6. Dependence of energies of complex solutions for XXZ spin chain and number of sites 
N = 15 on coupling constant A, which corresponds to cos <p. 

FIG. 7. Left hand sides of equations ( |4.2[) and (|4.3[) as a functions of k\ — &2 for N = 9 and 
N = 10. 
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